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CHAPTER I 


INTRODUCTION 

The desire for lighter and stronger structural materi- 
als has led to a rapidly growing use of fibrous composite 
materials in many modern applications . By combining the 
properties of more than one material having different 
strengths and weights, a structural member that is lighter 
and stronger than the isotropic material it replaces can be 
produced. An added benefit often gained is that due to this 
multiphase construction composites are frequently more 
insensitive to damage and damage growth than isotropic 
materials . 

Many different types of composite materials can be 
found in use today such as short fiber composites, woven 
cloth composites, and composites with randomly placed 
fibers. The list of matrix materials used to support these 
fibers is at least as long as the list of different fibers. 
In advanced structures such as airframes, however, most 
structural members are made of multiple layers, each con- 
sisting of continuous, parallel fibers embedded in a matrix. 
These layers, each of which is referred to as a lamina, are 
bonded together to form the structural member, typically 
called a laminate. 
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By orienting these continuous fibers in specific direc- 
tions the designer can adjust the material properties of the 
laminate to suit the loading to be encountered in use. Some 
composites having the predominate loading in a given direc- 
tion, such as landing struts and wing panels, are often con- 
structed of all unidirectional composite laminates. Others, 
which encounter a more randomly oriented load spectrum, are 
made by constructing a laminate from a series of properly 
oriented unidirectional lamina. 

In addition to the fiber orientation, advanced compos- 
ite materials are further grouped into two broad categories 
based on the type of matrix material used, these being plas- 
tic matrix composites and metal matrix composites. Unidi- 
rectional composites which have epoxy (plastic) matrices 
characteristically form large splits in the matrix between 
fibers near a defect such as a notch. Ductile metal matrix 
composites, on the other hand, exhibit matrix yielding 
between fibers near a notch instead of this splitting. It 
is found to be possible to have stable notch extension and 
broken fibers in composites which use ductile metals such as 
aluminum for the matrix material, presumably because of this 
large scale yielding. This is one question to be considered 
in the present study. 

For both unidirectional and multidirectional composites 
of the type discussed above, the fundamental component is 
the unidirectional lamina. It is thus essential that the 
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behavior of the lamina be clearly understood before the 
fracture behavior of the multiply laminate can be described. 
Many of the mathematical models which have been developed 
consider this fundamental lamina and the present experimen- 
tal study will attempt to relate the predictions of some of 
these models to actual laboratory measurements. 

Considerable analytical work has been done in modeling 
the unidirectional lamina with many of the recent mathemati- 
cal models developed at Clemson University by the authors 
[1,2]. These models, based on earlier work by Hedgepeth [3] 
in 1961, assume that the parallel fibers in the lamina carry 
all of the axial load and the matrix material transmits load 
to adjacent fibers by shear. An influence function techni- 
que was used by Hedgepeth to solve the lamina equilibrium 
equations. In later work, Hedgepeth and Van Dyke [4] were 
able to include matrix yielding and splitting in their 
model, although they were restricted by the solution techni- 
que to only one broken fiber in the lamina. Later models, 
such as the model of Zweben [5], extended Hedgepeth's model 
to include an arbitrary number of broken fibers while still 
allowing for matrix damage. Zweben, however, simplified the 
model considerably by grouping the fibers in three regions, 
one including all of the broken fibers, one having unbroken 
fibers which are influenced by the notch, and one group of 
unbroken fibers which are not influenced by the defect. 
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A limited amount of experimental results was available 
at the time of publication of [1], with perhaps the most 
complete being that of Awerbuch and Hahn [6] using narrow 
unidirectional boron/aluminum panels. These center-notched 
panels were loaded to failure so that their notched 
strengths could be determined. Extensive monitoring of the 
crack opening displacement (COD) was performed during the 
tests, along with some observations concerning the stable 
notch extension and matrix damage. This laminate damage was 
monitored through the use of radiographs and photomicro- 
graphs. The COD, defined as twice the vertical displacement 
of the center fiber, was monitored through the use of an 
interferometric displacement gage. 

Awerbuch and Hahn were able to predict the load-COD 
curves using a simplified analytical model which could 
account for some crack-tip yielding in the matrix. Results 
for a narrow range of initial notch widths were presented 
along with some interpretation of the relation between the 
observed failure modes and deformation characteristics. All 
of the tests performed by Awerbuch and Hahn, however, were 
carried out using narrow, 25.4 mm wide specimens. A finite 
width correction factor was used in the study to allow the 
predicted stresses for the infinite panel to be compared 
with the narrow specimens. These investigators assumed that 
the finite width correction factors available for isotropic 
plates would be applicable to their fibrous composite 
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specimens . But even with the finite width correction fac- 
tors , these investigators were limited to notch widths of 
12.7 mm (roughly 71 broken fibers) or less because of the 
narrow specimen width. All of the specimens used in the 
experiments were eight-ply; thus nothing was said concerning 
the effect of laminate thickness on the damage growth and 
failure mechanisms . 

Grandemange and Street [7] performed tests on double 
edge-notched boron/aluminum panels in an early attempt to 
investigate fracture mechanisms in fibrous composite materi- 
als. Their investigation included unidirectional laminates 
having a reasonably wide range of thicknesses, but a very 
narrow range on notch widths. The suggestion was made by 
these authors that the thickness of a specimen did affect 
the notched strength of an edge-notched panel. However, the 
strength of a panel having a given notch width varied only 
10 percent over the entire range of thicknesses tested. It 
is clear that a much wider range of notch widths must be 
tested before the thickness effect can be thoroughly 
explained. 

A detailed study of fracture of boron/ aluminum 
laminates containing unidirectional and + 45° fibers 
is presented by Poe and Sova in I8J . Reedy {9J 
considered thick (28 ply) unidirectional laminates and 
suggested in his paper that the thick specimens he used 
for his experiments should have been more notch 
sensitive than monolayer panels which he modeled 
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analytically. Yet in the results published by Reedy, the 
monolayer model predictions appeared to agree with the expe- 
rimental results for the thick specimens reasonably well. 
However, as in the work of Awerbuch and Hahn [6], only one 
thickness was used in the experiments, and the question of 
thickness effect on fracture behavior was not addressed. 

Reedy also noted the large amounts of longitudinal 

r 

matrix yielding predicted by the analytical model, and 
attempted to measure this yielding in the test laminate at 
various load levels using strain gage rosettes mounted at 
various points on the panel. The yield region predicted 
analytically appeared to agree fairly well with the strains 
measured in the laboratory, but since the strain gages used 
were 3 mm in width and thus covered nearly 20 fibers, there 
is some question as to the accuracy of these readings. The 
gages were clearly connected to some fibers that were broken 
and some that were not, thus the output of the gages needs a 
more careful analysis than is available from standard 
stress-strain relations. 

Other attempts have been made to describe the damage 
growth in fibrous composite panels during fracture. Rollins 
[ICj and Grenis and Levitt [11 ] attempted to relate the fiber 
breaks during a loading cycle to the acoustic emissions 
detected by transducers mounted on the specimens. It was 
found that the acoustic emission associated with a fiber 
break could be detected since these emissions were quite 
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energetic compared to those caused by matrix yielding and 
splitting. The points at which the various damage modes 
would occur in the loading sequence and the damage necessary 
to break a fiber were, however, never adequately explained. 

The need for additional experimental results became 
apparent during the final stages of the analytical modeling 
by Dharani, Jones, and Goree [ 12 ]. It was felt that a clear 
and logical description of the damage growth and fracture 
behavior in unidirectional laminates was provided by the 
micromechanical model. No detailed experimental studies 
describing this damage growth and fracture behavior were 
available, however, and much of the motivation for the work 
to be presented in this study was the need to observe very 
carefully the damage growth in a unidirectional composite 
and to compare the results with those predicted by the anal- 
ysis from [1] and [2]. 

The analysis given in [1] and [2] is valid for both 
brittle (plastic) matrix and ductile (metal) matrix compos- 
ites with the fracture behavior of a brittle matrix laminate 
such as graphite/epoxy being relatively simple to describe. 
At very low loads, the brittle matrix splits at the ends of 
the notch with the splits soon reaching the ends of the 
specimen for little increase in load. The specimen then 
fails in net-section tension with essentially no stress con- 
centration effect other than the loss in area due to the 
notch. 


7 



On the other hand the ductile matrix composite such as 
boron/aluminum exhibits large longitudinal yielding and 
transverse broken fibers and is thus much more complicated 
and difficult to model accurately. It was thus proposed 
that an experimental investigation be conducted using unidi- 
rectional boron/aluminum composite panels. It was known 
from past experience that finished boron/aluminum panels 
were available that possessed very uniform fiber size and 
spacing, and that the high quality of these finished speci- 
mens would make the tests much more valid. Very detailed 
records would be kept concerning the amount of matrix yield- 
ing or splitting as well as the number of broken fibers 
present in the laminate at various points in the loading. 

The results of these tests could then be compared to the 
damage zones predicted by the mathematical model. 

Unlike earlier tests, this proposed investigation would 
incorporate specimens of various widths and thicknesses 
(numbers of plies). Through the use of wider specimens, a 
much broader range of initial notch widths could be investi- 
gated than that of Awerbuch and Hahn [6]. These wider 
notchs could be cut in much wider specimens, reducing the 
need for the use of the finite width correction factor used 
in [6]. If the notchs in the proposed tests grew such that 
a finite width correction factor was needed, a much more 
accurate solution developed by Dharani [13 ] for the unidi- 
rectional composite laminate could be used. Also, the ques- 


8 



tion of thickness effect on damage growth proposed in [7] 
and [9] could conceivably be answered through the use of 
specimens having numbers of plies ranging from one to eight. 

As in some earlier investigations, the number of broken 
fibers present in the laminate at various load levels could 
be determined using X-rays to penetrate the specimens. The 
X-ray records would be obtained using a portable X-ray 
source which would allow radiography while the specimen is 
in the testing machine under load, rather than unloading and 
placing the specimen in an X-ray cabinet as was done in [6]. 
It was hoped that photoelastic coatings could be used to 
monitor the amount of matrix yielding present in the lami- 
nates, since this method had been used by others such as 
Yeow, Morris, and Brinson [14]. But since the thinnest 
available coatings were thicker than many of the proposed 
specimens, it was clear after some investigation that the 
photoelastic coating would not be able to detect the yield- 
ing between fibers with the desired accuracy. 

It was thus decided that a brittle lacquer surface 
coating would be the best method by which to measure the 
yielding present in the matrix between the fibers. It had 
already been shown by Chaturvedi and Agarwal [15] that a 
brittle coating would crack at a given threshold strain in a 
direction perpendicular to the maximum principal strain even 
on a unidirectional composite material. Because the brittle 
coating can be applied in a very thin layer it is sensitive 
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to local effects which occur over dimensions on the order of 


the fiber spacing. Thus the location and orientation of the 
cracks on the specimens in the proposed tests should enable 
the investigators to determine the extent of yielding pres- 
ent at selected load levels. 

One aspect of damage growth in composites that has not 
been investigated thoroughly is the effect of different 
matrix properties on the damage growth and fracture behav- 
ior. Only some very recent work by Reedy [16 ,17 ] has 
attempted to determine the result of using a stronger or 
weaker aluminum matrix material to hold the fibers. This 
work was again limited in scope, using only one notch width 
and narrow, 25.4 mm wide specimens. It was thus suggested 
that some of the proposed specimens have an aluminum matrix 
that had been heat-treated to a T6 condition instead of the 
fully annealed 6061 aluminum matrix in the remainder of the 
panels. This stronger matrix should alter the fracture 
behavior, since it is through the matrix that the load is 
transferred from the broken fibers to the remaining unbroken 
ones. Furthermore, unlike the work in [!]>•■] and [17]/ a wide 
range of notch widths and panel thicknesses could be tested 
to provide much more information about the effect of matrix 
properties . 

A series of tests using a wide range of specimens hav- 
ing different notch widths, overall widths, numbers of 
plies, and matrix properties was thus proposed. A complete 
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description of the materials and test procedure used in the 
tests is described in the next section. 

Certain commercial materials are identified in this 
paper in order to specify adequately which materials were 
investigated in the research effort. In no case does such 
identification imply recommendation or endorsement of the 
product by NASA, nor does it imply that the materials are 
necessarily the only ones or the best ones available for the 
purpose . 
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CHAPTER II 


EXPERIMENTAL PROCEDURE 
Materials 

The specimens used in this investigation were all uni- 
directional boron/aluminum composite panels of various 
widths and thicknesses. The panels, manufactured by Amer- 
com, Inc. in Chatsworth, California, consisted of 0.142 mm 
(5.6 mil) diameter boron fibers in a 6061 aluminum matrix. 
Some of the panels were heat-treated after the initial manu- 
facture such that the aluminum matrix would possess the 
properties of standard 6061-T6 aluminum instead of an ann- 
ealed 6061-0 aluminum matrix. Equal numbers of one-, two-, 
four-, and eight-ply specimens were fabricated, with a 
thickness of approximately 0.178 mm per ply. A complete 
listing of the test specimens used in this study can be 
found in Table (I). It should be noted that a set of speci- 
mens is defined as one single-ply specimen, one two-ply 
specimen, one four-ply specimen, and one eight-ply specimen, 
all having the same notch size and laminate width. Slightly 
thicker aluminum foil was used on the single-ply panels in 
order to increase out-of-plane stiffness. Thus the fiber 
volume fraction for the single-ply panels was around 35 per- 
cent, while the volume fraction for all other specimens was 
very close to 50 percent. 
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Table I. List of specimens used in experimental 
investigation . 


Panel Width (W) 

Notch Width 

(2a) 


* 

Matrix Temper 


25 . 4 mm 

1.59 

mm 



1 set 6061-0 


25 . 4 mm 

3.18 

mm 

1 

set 

6061-0, 1 set 

6061-T6 

25 . 4 mm 

4.76 

mm 

1 

set 

6061-0, 1 set 

6061-T6 

25 . 4 mm 

6.35 

mm 



1 set 6061-0 


50.8 mm 

1.59 

mm 



1 set 6061-0 


50.8 mm 

6.35 

mm 

1 

set 

6061-0, 1 set 

6061-T6 

50.8 mm 

9.52 

mm 

1 

set 

6061-0, 1 set 

6061-T6 

50 . 8 mm 

15.9 

mm 



1 set 6061-0 


73.1 mm 

3.18 

mm 



1 set 6061-0 


73.1 mm 

12.7 

mm 

1 

set 

6061-0, 1 set 

6061-T6 

73.1 mm 

15.9 

mm 

1 

set 

6061-0, 1 set 

6061-T6 

73.1 mm 

19.1 

mm 



1 set 6061-0 


101.6 mm 

15.9 

mm 



1 set 6061-0 


101.6 mm 

19.1 

mm 

1 

set 

6061-0, 1 set 

6061-T6 

101.6 mm 

22.2 

mm 



1 set 6061-0 


101.6 mm 

25.4 

mm 

1 

set 

6061-0, 1 set 

6061-T6 


One set equals one single-ply panel, one two-ply 
panel, one four-ply panel, and one eight-ply panel. 
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Specimens having widths of 25.4, 50 . 8 , . is . 1 f and 100.6 
mm were cut from large panels obtained from the manufac- 
turer. These cuts were parallel to the fibers and were made 
with a diamond saw. The edges were later ground to ensure 
minimum surface damage due to cutting. All of these experi- 
mental specimens were approximately 30 cm in length. Center 
notches of various widths were then introduced by E.D.M. 
(Electrostatic Discharge Machining). Normally the width of 
the notch did not exceed one-quarter of the overall width of 
the panel. 

Finally, the specimens were carefully cleaned and pol- 
ished to eliminate further the possibility of surface flaws 
that might affect the testing. Two strain gages were 
mounted on each specimen at points well below the notch so 
that the remote strain could be measured and recorded. The 
strain gages also enabled the experimenters to determine if 
the panel was indeed being subjected to a uniform remote 
stress. The last step in the preparation was to epoxy alu- 
minum doublers to each end of each panel so that the testing 
machine grips would not damage the ends of the specimens. A 
representation of a typical test specimen is shown in Figure 
( 1 ). 
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Figure 1 . 


Typical Test Specimen. 
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Radiographic Technique 


The radiographs used in these experiments were taken by 
exposing the specimens to low level X-rays for a specified 
amount of time. The source of these X-rays was a Model MTK 
140Be X-ray machine manufactured by the Philips Company of 
West Germany. This is a portable unit and has the signifi- 
cant advantage that the radiographs can be taken while the 
specimens are in the testing machine and under load. 

Much time and experimentation was involved in determin- 
ing the best exposure time and film to focus distance (FFD) 
for each specimen thickness. The X-ray source used was 
capable of varying both the voltage and the current which 
produce the X-rays. It was quickly discovered that each of 
the variable parameters must be chosen correctly if an 
acceptable degree of resolution was to be realized. Many 
test shots were taken of a given laminate at various cur- 
rents, voltages, distances, and exposure times. The nega- 
tives from these tests were then used to make radiographic 
prints to determine the correct exposure parameters. 

It was found that the single-ply and the two-ply speci- 
mens could be x-rayed at an FFD of 50 cm with a current of 5 
mA. On the thicker four-ply and eight-ply specimens, the 
best radiographs were taken at an FFD of 80 cm with a cur- 
rent of 5 mA. The higher FFD serves to place the X-ray 
focal point further away from the object, thus making the 
X-ray seem more like a point source; that is, the rays are 


16 



almost parallel when they strike the specimen. The voltages 
and exposure times used varied depending on the thickness of 
the panel. A graph of current-exposure time product versus 
the specimen thickness for various voltages is shown in Fig- 
ure (2). Any X-ray exposure taken with data from this graph 
was found to give excellent resolution on the radiographic 
print. As an example, the four-ply panels, with a thickness 
of 0.712 mm, were x-rayed for three minutes at 30 kV with a 
tube current of 5 mA, giving the required product of fifteen 
mA-min as determined from the graph. 

The image produced from the X-ray is captured by simply 
securing a sheet of Type 55 Polaroid film to the back of the 
specimen at various load levels during the test. It must be 
emphasized that the film is not placed in a film holder, 
such as a Polaroid Model 455, which has an intensifying 
screen. The intensifying screen is not needed for this type 
of exposure; in fact, use of the film holder actually made 
the resolution worse on the radiograph. Thus the sheet of 
film is taped to the specimen encased only in the light 
paper provided by the manufacturer. 

The number of X-rays absorbed or transmitted by a mate- 
rial is dependent on its atomic number. For example, an 
X-ray of a graphite/epoxy panel will produce no distinct 
fiber pattern on a radiograph since both the fibers and the 
matrix are composed of essentially carbon. In fact, boron 
(atomic number =5) and aluminum (atomic number =13) 
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CURRENT-EXPOSURE TIME PRODUCT (mA-min) 



Figure 2. X-ray Exposure Chart. 
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transmit amounts of X-rays that are so close that it is dif- 
ficult to distinguish between the two of them on a radio- 
graph. The boron fibers commonly available, however, are 
manufactured by depositing the boron on a thin tungsten 
filament, and this tungsten (atomic number = 74) core will 
appear as a very dark black line on a radiograph. Thus the 
tungsten core of each fiber can be easily identified on the 
radiographic print. 

A typical radiograph of a single-ply boron/aluminum 
panel is shown in Figure (3). The dark lines are the tung- 
sten cores of the boron fibers, while the light area around 
these cores is the boron itself. The medium gray area 
between the fibers is the aluminum matrix. Thus as the 
specimen is loaded and fibers begin to break, it is a simple 
matter to count the number of broken fibers ahead of the 
notch at various load levels. It is also possible to obtain 
a very accurate measure of the crack opening displacement 
(COD) by merely scaling the appropriate fiber displacements 
directly from the radiographic print made in the darkroom. 
COD measurements were also taken with a clip gage which is 
specifically designed to measure crack opening displacement 
for notched specimens. The output of this gage can be care- 
fully calibrated such that the strain output from the gage 
can be converted to a displacement in millimeters. 
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Figure 3. Radiograph of single-ply boron/aluminum composite panel 
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Brittle Coating Technique 

A brittle lacquer was applied to each specimen to meas- 
ure the amount of longitudinal yielding between the last 
broken fiber and the first unbroken fiber at the end of the 
notch. The brittle lacquer used in these tests was the com- 
mercially available Tenslac, manufactured by the Measure- 
ments Group in Raleigh, North Carolina. This lacquer will 
crack when a certain threshold strain is reached in the 
boron/aluminum panel. This threshold strain can be affected 
by the coating thickness as well as the temperature and 
humidity of the test environment. 

As previously mentioned, each specimen was carefully 
cleaned and polished to provide a specimen surface that was 
free from surface flaws. The polished surface also provided 
a highly reflective surface below the brittle lacquer. This 
made the cracks reflect the light much more clearly so that 
the extent of cracks could be observed and photographed. It 
was noted that the ability to see a given crack in the coat- 
ing was strongly affected by the angle at which the light 
approached the surface of the specimen. The cracks could 
not be seen if the light approached the surface along a line 
normal to that surface. Rather, the best results were 
obtained when the approaching light made an angle of approx- 
imately 30 degrees with the surface normal. 

The brittle lacquer was applied to each specimen in up 
to twenty very thin coats to minimize bubbling and provide 
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the desired "orange-peel" finish. The final thickness of 
the lacquer coating was approximately 0.075 millimeters 
which gave a threshold strain of 1200 to 1800 microstrain 
depending on test conditions. The threshold strain for each 
specimen tested was measured by using strain gages that were 
mounted on each panel on the opposite side from the coating. 
The remote strain at which the cracks appeared directly 
opposite the gages was taken as the threshold strain for 
that specimen. This procedure eliminated the need of pre- 
paring calibration bars or even assuming a nominal threshold 
strain. Rather, the exact threshold strain at the time of 
the test for each specimen could be accurately measured and 
recorded. The brittle lacquer coating also gave a very 
clear picture of the stress state at the remote sections of 
the specimens. If the remote stress was truly uniform, the 
cracks in the coating would be straight and perpendicular to 
the fibers outside the region affected by the defect. The 
output from the strain gages was measured through the use of 
a Vishay VE-20 Strain Indicator and a Vishay VE-21 Switching 
and Balancing Unit. 
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General Testing Procedure 


r 

After allowing the brittle coating to dry for 18 to 24 
hours, the specimens were loaded statically with a Baldwin 
tensile testing machine. A photograph of the test set-up is 
shown in Figure (4). A baseline X-ray was taken of each 
panel before the loading cycle was begun so that any damage 
caused by machining or polishing could be noted. Then as 
the loading sequence began, photographs were taken of the 
cracks in the coating at various load intervals with a Nikon 
FM camera equipped with a close-up lens. A bright light was 
used to illuminate the crack pattern so that a small lens 
aperture could be set to increase depth-of-f ield. . 

X-rays were taken in the manner described earlier at 
selected load levels so that the number of broken fibers in 
the specimen could be determined. The majority of the X-ray 
exposures were taken in the last half of the load cycle 
since the first fiber breaks occurred at approximately 50 
percent of the failure load. The crack opening displacement 
was measured both from the radiographs and with a clip gage 
specifically designed to measure COD. 

A total of ten to twelve X-rays and twelve to fifteen 
brittle coating photographs were made during each test. The 
strain output from each gage was recorded at numerous load 
levels as were the load readings from the Baldwin machine. 
This testing procedure was repeated for boron/aluminum pan- 
els of many different widths, thicknesses, and initial notch 
widths . 
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Figure 4 


Experimental Test Set-up 


CHAPTER III 


ANALYSIS 

Mathema tical Model Description 
Though the main emphasis of this paper is the experi- 
mental investigation, it is important that the mathematical 
model which motivated the present study and the assumptions 
made therein are clearly understood. Thus at this point, a 
brief description of the equations involved in the deriva- 
tion of the shear- lag model will be presented. Following 
Goree and Gross [1], the laminate is modeled as a two-dimen- 
sional region, having a single row of parallel, identical, 
equally spaced fibers, separated by matrix as shown in Fig- 
ure (5). The initial damage is taken to consist of an arbi- 
trary number of broken fibers such that all breaks lie along 
the x-axis. The longitudinal matrix damage is introduced at 
the end of the initial notch as in [1]. The additional 
notch tip transverse damage consists of an arbitrary number 
of broken fibers which are constrained by the adjoining 
matrix and/or the unbroken fibers through the thickness. 
These fibers in the transverse damage zone will be referred 
to as constrained fibers. It is mathematically untractable 
to account precisely for the distribution of fiber damage 
occurring in the region ahead of the notch. The model 
assumes that all the breaks occur on the x-axis and accounts 
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Figure 5. Two-dimensional unidirectional lamina with broken 

fibers, longitudinal matrix splitting and yielding, 
and transverse fiber and matrix damage (first 
quadrant) . 
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for the reduction in stiffness of the transverse 

damage 20ne by defining a stiffness coefficient, X, such 

that 


- stress in the constrained fiber 
® ~ stress in the first unbroken fiber * 

The free-body diagram shown in Figure (6) is considered 
with additional conditions for the last broken fiber, 
denoted by n = N, and y < L, to account for the longitudinal 
damage taken as 

Ti N+l = ' T o <Y " l> ' (1) 


where 


<y - £> = 1, y > £, and 

<y - £> = 0, y < £ . (2 ) 

L is defined as the total damaged length (see Figure 5), £ 

as the split length, and t as the matrix yield stress. 

o 

Yielding is assumed to occur when the matrix reaches the 
yield strain, X Q . Splitting occurs at a multiple of X Q as 
given by the particular matrix material. 

The equilibrium equations in the longitudinal and 
transverse directions for all fibers n, with the exception 
of N and N + 1 when y < L, are then 


27 






0 


(3) 


A F dg F>n _ . . . = 

t dy *n+l 'n 


and 


°M^n+l 


M 1 n 


+ 


h d , , 

{t| n+l 


+ t| } = 0. 
'n 


(4) 


For fiber N, when y < L, t I n+ ^ = - t <y - Z> , and the 
equilibrium equations are 


^F da F | N 
t dy 


t <y 
o 


Z> 



(5) 


and 


°M 1 N+l °m'N ' 2 dy 


'n + <y " e> " "N> = < 6 > 


For fiber N+l, when y < L, T I n+ 1 = 
the equilibrium equations are 


t q <y - Z > , and 


tZ do F | N+l + . + , <y - £> = 0, 

t dy 'n+1 o y 


(V) 


and 


°M*N+2 " °M I N+1 + I'dy { ' c| N+2 “ T o <y ‘ a>} °* 

The stress-displacement relations common to the shear- 
lag assumption. 
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( 9 ) 


dv 

°F 'n = E F*^T' 
T, n+1 = G M (v n+l 
°M ^ n+1 = E M (u n+l 


- ' 
- u n )/h 


and 


( 10 ) 

( 11 ) 


are introduced into the equilibrium equations, leading to 
the replacement of the three stress quantities by the axial 
and transverse displacements. These assumptions, however, 
simplify the equilibrium equations by removing the tran- 
sverse displacement dependence from the longitudinal equi- 
librium equation. The fiber stress and the matrix shear 
stress can then be determined without solving the trans^ 
verse equilibrium equation. Only the longitudinal equilib- 
rium equation will now be considered. For all fibers, 
except N and N+1 where y < L, the longitudinal equilibrium 
equation becomes 


E F A F h d v n 
G M t dy 2 


+ V 


n+1 


2v + v , = 0. 
n n-1 


( 12 ) 


For fiber N where y < L, 


E^Ar-h d v N h 


G-wt 

M 


(13) 


For fiber N+1 where y < L, 


Mf 11 d V N+1 
G M t dy 2 


+ V 


N+2 


v „ + — x <v 
V N+1 G m o Y 


- Z> = 0 . 


(14) 
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It is convenient to use the non-dimensional forms of 
Equations (12), (13), and (14). Noting the coefficient of 

the second derivative term in the above equations, define 


y = 


E F A p h 




" and °Fln = °eo°n 


The other required quantities, such as the axial displace- 
ment (v ) and the matrix yield stress (t ), can also be 
n o 

written in their non-dimensional forms. From the defini- 
tions of y and Opl n it follows that 



A v h 

1% t h 


V = 0 

r 

> 

3 f 

II 

> 

wher 

n oo 

E p G M t 

n T^ M n 



g m a f 

% T 

T or a = 

Epht’ 

T — O 

o 00 

Epht 

°° T 

g m a f_ 


% T 


T. 


dV 


I = a 
F 'n o 


n 


dn 


o 

T 


Epht 


g m a f 


% dV. 


n 


dn ' 


IVfU „ . _^o . _ .. . 

T n 0 °°|E F htJ * V n V n-1^ ~ V n ~ V n-1^' 


L = 


|E F A F h| % 


and ft = 


Ep A ph 

G M C 


% 


3- 


M15) 


The quantities n, t, o n , V , o, and 3 are non-dimensional, 
whereas Ep, Ap, t, L, and d represent the actual fiber modu- 
lus, fiber area, lamina thickness, and damage zone dimen- 
sions, respectively. 
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The resulting equilibrium equations in non-dimensional 
form are given by 


d 2 V 

—r + V n+1 - 2v n + Vl = °' 


dn' 

d 2 V. 


N 


dn 

d 2 v. 


f + V N + 1 - 2V N + V N-1 = - f(T »' and 


dn 


T 1 * V N + 2 - 2V N + 1 + V N = f( ">- 


(16) 


(17) 


(18) 


where f(Tj) is a new unknown function such that 


f(n) = v n “ V N+1 ~ T 0 <T * " e>/ ^ < a, and 

f (n) =0, T1 > a. (19) 


These differential-difference equations may be reduced to 
differential equations by introducing the even valued trans- 
form as 

oo 

V(T,,6) = V_(n)/2 + 1 v_(n) cos (n0 ) , (20) 

0 n=l n 

from which 

2 * - 

V (u) =- / V(n,0) cos(n0)d0 . (21) 

n k o 

Making use of the above transformation and the orthogo- 
nality property of the circular functions, the three 
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equilibrium equations may be written as one equation valid 
for all values of n and 11 as 


TT 2— 

— / { ^ -2 [ 1 - cos(6)]V} cos(n0)d8 

71 0 dri 


2 

IT 


X 


TT 

<a-ii> J f (ti ) {cos [ (N+l ) B ] 

0 


cos(N0)} cos(n8)d0. 

( 22 ) 


This equation is of the form 


2 17 

— J F(ti, 0) cos(n0)d0 = 0 for all ti and n. 

77 0 


Noting the definition of V(ti, 0) in Equations (20) and 
(21) it is seen that the function F(n,0) is even valued in 0 
and therefore, if the integrand is to vanish for all n, the 
function F(n,0) must be zero. The single equation specify- 
ing V(n,0) is then 



e 2 v 


<0C-T| > D f ( Tl ) 


(23) 


where 


2 2 

6 = 2[1 - cos(0)] = 4 sin (0/2), and 

2 

D = cos(N0) - cos[(N+l)0]. 
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The solution to the problem of vanishing stresses and 


displacements at infinity and uniform compression on the 
ends of the broken fibers will now be sought. The complete 
solution is obtained by adding the results corresponding to 
uniform axial stress and no broken fibers to the following 
solution. The appropriate boundary conditions are thus 


V n =0 as ti (24) 

dv n 

=0 as ti -*• and (25) 


V n = 0 for ti = 0, for unbroken fibers, (26) 


dV n — 

-g— = o n = -1, for broken fibers ( n=0 to n=N), (27) 


= -1 + Z 


, for constrained fibers 


(n=N+l to n=M) 


The complete solution to Equation (23) satisfying van- 
ishing stresses and displacements at infinity is given by 


V(n,9) = A( 8 ) e" 5n + / (sinh[ 8 (n-t) ] <o-t> f(t)dt, 

n 

(29) 


where the unknown functions are A(0) and f(t). The remain- 
ing boundary conditions give 
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dV ( 0 ) 9 7T 9 a 

— = ± j {-6A(0) + D J cosh(6t)f (t)dt}cos(n9)d6 = -1, 


n = 1, 2, . . . , N, 


(30) 


dV (0) o 71 o a 

— 5 =- I { - 6A ( 0 ) + D / cosh(6t)f (t)dt) X 

ar ) " n n 


cos(n0)d9 = -1 + n = N+l, . .., M, 


(31) 


2 Tf D 2 a 

V (0) = — / {A( 0 ) — 7- I sinh( 6t ) f ( t ) dt} cos(n0)d0 = 0, 
n 7T o 0 0 

n = M+l, . . . , (32) 


Equation (32) is solved exactly by taking 


2 (X oo 

A(0) - f sinh(6t)f ( t)dt = I B_ cos(m0) 
6 0 m=0 m 


(33) 


where m is a broken fiber index and the 13 are constants. 

m 

The stress in the first unbroken fiber, a F I M+l ‘*' s g ^ ven ky 

°^M+1 2 71 2 a 

- ; - = ■£ I {-6A( 9 ) + D J cosh(5t)f (t)dt) X 
ar) TT n n 


cos[ (M+l )0 ]d8 = -1 + ° F ! M+l * 


(34) 


Using Equations (33) and (34) in Equations (30) and 
(31), A( 0 ) and ma ^ e li m i nate d, resulting in M+l 

algebraic equations for M+l constants B m in terms of f(n) 
which is, as yet, unknown. For longitudinal matrix damage. 
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Equations (30) and (31) must be supplemented by the condi 
tions that 

f(n) = g(n) - x <i) - P>, ti < a, and (35) 

g(ct) = T , (36) 

where 


g(Tl) ” V N V N+1 * 


From Equations (28) and (33), A(6) may be eliminated to 
obtain V(ti, 8) in terms of constants B m and the unknown func- 
tion f(t). Recalling the relation between V(ti, 0) and V n (n) / 
an expression can be obtained for the axial fiber displace- 
ment V (ti) as 
n v 


0 7r M 

V (ti) = — / e -01 ' I B cos(m0)cos(n8)d0 
n u 0 m=0 m 

1 a 

+ 4 J f(t) { C ( ] t-T) | ) - C ( t+Tl ) > dt, 

2 o n n 


where 


n 


(£) 


_ 2 / 

7T * 5 




cos(n8)d0. 


(37) 
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Equations (30) and (31) then become 


9 7T M 0 a _*«- 

— / {-6 Z B cos(m0) + D z J e” c g(t)dt 
77 0 m=0 m 0 


- D 2 x / e < ^ t dt} cos(n8)d0 = -1, 


n = 0, 1, . . . , N, 


( 38 ) 


and 


9 7T M n a 

- / {-6 Z B_ cos(m0) + D / e g(t)dt 

77 0 m=0 m 0 


ct 

- D 2 7 / e _<St dt]{cos(n0 ) - Ycos [ (M+l ) 0 ] >d0 = -1 + K, 

3 

n = N+l, N+2, . . . , M. (39) 


Equations (36) and (37) give 


? it M 

g(ri) =— / e Z B_ cos(m0) (cos(n0) 

77 0 m=0 m 


- cos [ (N+l ) 0 ] ) d0 + 


+ ^ /“g(t) {C N ( | t-n | ) - c N (t+n) - c N+1 ( | t-n | ) 


+ c N+ 1 (t+n)} dt 


— a 

T 


j S {c N (|t-n|) - c^j (t+n) 

3 


- ^"N+l ( ^ ^-~ T i I ) 


+ C N+1 (t+n)> dt. 


( 40 ) 
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which is a Fredholm integral equation of the second kind. 

The last condition that must be satisfied is that 

g(«) =~ (4i) 

The series equation and the integral equation are cou- 
pled by the Fourier coefficients and the function g ( n ) and 
may be solved on a digital computer using a Gauss-Legendre 
quadrature scheme. Recent improvements in the solution 
scheme have reduced the computation time to less than two 
minutes, even for a very large number of broken fibers and 
large matrix damage. 

Brittle Coating Interpretation 

The mathematical model described in the previous sec- 
tion gives the fiber displacements and stresses throughout 
the laminate as well as the longitudinal yield region 
between the last broken fiber and the first unbroken fiber. 
An important step in relating the behavior predicted by the 
model to the measured laboratory values is thus to compare 
the measured longitudinal yield zone with the predicted 
value of the yield length from the model. The difficulty in 
this comparison is the interpretation of the crack pattern 
found in the brittle coating photographs. 

Using the fiber stresses and displacements from the 
mathematical model it is possible to approximate the surface 
stress state in the aluminum matrix between the fibers. 
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Since the thin coating of brittle lacquer is sensitive to 
strains in this small region, the angle of the cracks in the 
coating should correspond to the principal strain (or 
stress) angles calculated from a simple Mohr's circle analy- 
sis for the aluminum. It should be noted that the matrix 
normal stress in the vertical direction was omitted in the 
shear-lag analysis for the fiber stresses since this stress 
is small compared to the fiber stress. However, it is 
appropriate to now consider this matrix stress since it is 
part of the stress state in the matrix which influences the 
crack pattern in the coating. 

A good approximation for the average matrix normal 
strain in the vertical direction can be calculated by aver- 
aging the strains in the two adjacent fibers. 


Vn = 


a 

oo 

TK 


dV 


n 


dV. 


n+1 


dti 


dn 


(42) 


The horizontal matrix strain is assumed due to a Pois- 
son contraction alone, and is thus derived by multiplying 
Equation (42) by the Poisson ratio for aluminum, 0.3. The 
stresses in the matrix are then a = ®M e y an< d ° x = 0. 
Finally, the matrix shear stress is calculated from Equation 
(15) as 


T = — — 

n - 


< V n - V n + 1> 


(43) 
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Thus a typical element as shown in Figure (7) below has 
a known state of stress and the angle through which the ele- 


ment must rotate to reach the 
be calculated from 

2t xv 

tan (26 ) = - _ y g 

x y 



angle of principal stress can 



Figure 7. Mohr's Circle Analysis. 

The growth of longitudinal damage predicted from the 
above analysis can now be compared to the growth from actual 
specimen photographs. The longitudinal crack pattern for a 
specimen having an overall width of 73.2 mm and an initial 
notch width of 12.7 mm will be examined at various load lev- 
els and compared to predicted strain data from the model. 

The first load step examined is very early in the cycle 
before any longitudinal yielding has begun. At an axial 
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load of 800 N the model predicted that very few cracks 
should form, since the strain falls below the threshold 
strain of 1500 ye at a y-value of approximately 1 mm. The 
model further predicts that these first cracks that form in 
the coating will form at an angle of 37 degrees from the 
horizontal. Table (II) shows the predicted data from the 
model at the load levels to be examined in this analysis. 
Figure (8) shows the appearance of the coating at a load of 
800 N. It is clear that these first cracks, forming an 
angle measured to be 32 degrees, appear directly above the 
end of the notch between the last broken fiber and the first 
unbroken fiber. Moreover, the uppermost crack forms at a 
y-value on the order of vthat predicted by the strain values 
computed from the analysis. 

At a load of 3000 N the model predicts that approxi- 
mately 2.5 mm of longitudinal matrix will yield and that the 
principal strain in the matrix will exceed the brittle lac- 
quer threshold strain value for y-values less than 5 mm, 
with the last predicted crack forming at an angle of 30 
degrees. Examination of Figure (9) reveals that the upper- 
most crack, at a y-value of 6 mm, forms an angle of 29 
degrees with the horizontal. In this figure it can also be 
seen that a region of very intense cracking has formed, con- 
sisting of short, closely spaced cracks and lacquer crazing. 
The length of the region containing these short cracks cor- 
responds very closely to the yield length of 2.5 mm pre- 
dicted by the model. 


41 



The length of the longitudinal matrix yield region at 
4500 N is predicted by the model to be approximately 5 mm, 
and the uppermost crack in the coating at this load should 
form 8 mm above the notch at an angle of 29 degrees. Figure 
(10) shows that the uppermost crack is present 12 mm above 
the notch at an angle of 27 degrees and that the region of 
severe cracking and crazing has grown to approximately 5 mm. 

Finally, in Figure (11) the uppermost crack forms 14 mm 
above the notch at an angle of 20 degrees with the horizon- 
tal, and the region of severe cracking has grown to 10 mm at 
this load level of 6500 N. The predicted location and 
orientation of the last crack at this load level are 15 mm 
above the notch at an angle of 22 degrees with the horizon- 
tal. The yield length predicted by the mathematical model 
at this load level is 10 mm. Thus the stresses and strains 
in the matrix actually measured in the laboratory appear to 
agree reasonably well with those predicted by the mathemati- 
cal model. The growth of the crack region as well as the 
severely cracked region is shown in graphical form in Figure 
(12). It is clear that the brittle coating provides a very 
good qualitative picture as well as some quantitative data 
concerning the yield region at the end of the notch and its 
growth as the load is increased. 

It should be noted that the cracking discussed above 
occurs relatively early in the load cycle before the first 
fiber break has occurred. As fibers begin to break in the 
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laminate, the loading is stopped and restarted each time an 
X-ray exposure is made. As is the case with all brittle 
lacquers, the load must be removed completely each time the 
loading is stopped if the threshold is to remain valid. Of 
course, this would greatly alter the stiffness of the previ- 
ously yielded aluminum matrix and thus change the behavior 
of the laminate. The loading cycle was thus . continued after 
each X-ray exposure without first removing the load. 

The brittle coating, while providing actual numerical 
approximations in the early part of the load cycle before 
the first fiber has broken, can provide only qualitative 
results later in the cycle. It is felt, however, that the 
severely cracked region is closely associated with yielding 
in the matrix. Several tests were performed using pure alu- 
minum notched panels and in each case the region(s) in which 
yielding was known to occur exhibited this same intense 
cracking. The region of intense cracking and flaking is 
thus used as indicator of matrix yielding and will be refer- 
red to as the severe damage length in the remainder of this 
paper . 


43 



Table II. Brittle Coating Interpretation Data. 

Uppermost Crack Location Uppermost Crack Angle Severe Damage Region 


Load 

Predicted 

Observed 

Predicted 

Observed 

Predicted 

Observed 

800 N 

1 mm 

1 mm 

37° 

32° 

0 

0 

3000 N 

4 mm 

5 mm 

30° 

29° 

2 . 5 mm 

2. 5 mm 

4500 N 

8 mm 

11 mm 

29° 

27 ° 

5 mm 

5 mm 

6500 N 

14 mm 

18 mm 

22° 

20° 

10 mm 

10 mm 




r r an axial load 
rrack pattern for an a 
o Brittle coating crack P 
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Figure 10. Brittle coating crack pattern for an axial load 
of 4500 N. 
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Figure 11. Brittle coating crack pattern for an axial load 
of 6500 N. 
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CHAPTER IV 


RESULTS 

Damage Growth Sequence 

The main objective of this paper is to determine accu- 
rately the properties of boron/aluminum laminates and then 
compare the predicted values to those measured in the labo- 
ratory using actual specimens. The experimental tests 
should provide values for the matrix yield length, all fiber 
displacements, fiber stresses, and stable notch extension. 
Material properties and laminate dimensions from handbooks 
and manufacturer's specifications should allow the experi- 
menters to convert the non-dimensional predicted values from 
the model to real values suitable for comparison with expe- 
rimental results. 

As previously discussed, the use of radiographs and 
brittle coating photographs will present a very detailed 
picture of the growth of damage during a loading cycle of 
any unidirectional boron/aluminum laminate. It is therefore 
appropriate to begin the examination of the results of this 
study with a step-by-step presentation of a typical set of 
photographs of one of the experimental specimens. The spec- 
imen chosen for this analysis is typical of all of the test 
panels, exhibiting both matrix yielding and stable notch 
extension in the form of broken fibers ahead of the notch. 
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The set of photographs and radiographs show the pro- 
gression of damage in a two-ply boron/aluminum laminate 
which is 73.1 mm wide. A 19.05 mm wide notch (109 broken 
fibers per ply) was cut perpendicular to the fiber direction 
in the center of the panel before the testing procedure was 
begun. Brittle coating photographs and X-ray exposures were 
made at selected intervals in the loading sequence. The 
series of radiographs and photographs is presented in Fig- 
ures (13) through (28) beginning at the end of this section. 

As previously mentioned, a baseline X-ray was taken of 
each specimen before the loading cycle was begun to detect 
any damage which may have been caused by machining or pol- 
ishing. The baseline X-ray for the two-ply specimen is 
shown in Figure (13). The right end of the notch can be 
seen in this radiograph, along with the fiber cores which 
run perpendicular to the notch. The initial notch opening 
was 0.381 mm; thus crack opening displacement at each load 
level could be measured from the radiographs once the scale 
factor had been determined. It is clear that there are more 
fibers visible in this radiograph than in the earlier sin- 
gle-ply radiograph. This indicates that the fibers in each 
ply, as expected, are not perfectly aligned with the fibers 
in adjacent plies. It can be seen, however, that there are 
usually two fiber cores very close together on the radio- 
graph, indicating the two fibers in adjacent plies. 
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In the first few photographs, the cracks in the brittle 
coating grow in the longitudinal direction, even though no 
fiber breaks have occurred. In Figure (17) it can be seen 
that the first fiber breaks have occurred, indicated by the 
breaks in the fiber cores near the end of the notch. It 
should be noted that from this point there will be no longi- 
tudinal growth of the severe damage region, since the first 
previously unbroken fiber is now free to displace. The pro- 
gression of damage, indicated by severe cracking and flaking 
of the brittle lacquer, will now be in the transverse direc- 
tion as more fibers continue to break. 

First fiber failure typically occurred at 55-60 percent 
of the ultimate failure load for all laminates tested, inde- 
pendent of the laminate width, thickness, initial notch 
width, or the matrix temper (6061-0 or 6061-T6) . Subsequent 
fiber failure normally occurs at 75 to 80 percent ultimate, 
with visible notch extension occurring at roughly 95 percent 
ultimate. This behavior is described in Table (III). The 
present specimen shown in Figures (13) through (28) were 
chosen primarily because of the quality of the brittle coat- 
ing photographs and the radiographs, even though in this 
sequence the second fiber failure does not occur until 93 
percent of the ultimate laminate failure load. 

In Figure (24) the severe damage region in the coating 
has grown transversely from the notch, clearly indicating 
the existence of more than one longitudinal damage region in 
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Table III. Loads at which subsequent fiber breakage occurs for 73.1 min B/Al 
panels with 15«9 mm notches (89 broken fibers per ply). 


U1 

u> 


Matrix 



1- piy 

2 -ply 

4 -ply 

g-piy. 



r 

First fiber break 

55 % 

50 % 

60 % 

55 °fo 

6061 -0 


Second fiber break 

85 fo 

80 % 

83 % 

72 



Beginning of visible 

98 % 

97 % 

94 °]o 

95 



notch extension 







f First fiber break 

56 1o 

60 % 

59 # 

55 % 

6061 -T 6 

< 

Second fiber break 

83 % 

81 % 

79 °Io 

78 % 



Beginning of visible 

98 % 

94 % 

97 % 

93 % 


^ notch extension 


* 


These loads are given as a percentage of the ultimate notched failure 
load for each particular specimen,, 






the matrix. It is this phenomenon that allows the COD to 
grow at a faster rate than predicted by the model, since the 
model at the present time can include only one yield region. 
The last figures in the damage sequence indicate that the 
transverse damage region continues to grow until the failure 
load of the laminate is reached. It should be noted that 
the extension of the notch in the form of broken fibers does 
not occur in a straight line parallel to the notch. Rather, 
the crack grows in a jagged fashion with the fiber breaks 
occurring at different vertical locations as the crack grows 
and more fiber breaks are observed. The fibers in this 
region are thus broken but still able to carry some reduced 
amount of load since the matrix around them is still intact. 
This phenomenon exactly parallels the assumptions made for 
the constrained fibers in the previous section, indicating 
the physical basis for the constrained fiber approximation. 
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Figure 14. Brittle lacquer crack pattern at 2670 N showing 
formation of first cracks in the coating. 
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Figure 15. Radiograph of two-ply panel at 5340 N (35 t failure load). 


Figure 16. 


Brittle lacquer crack pattern at 5340 N. 
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Figure 17. Radiograph of two-ply panel at 8010 N (53% failure load) showing first 
fiber break occurring at the upper, right-hand portion of the notch. 



Figure 18. Brittle lacquer crack pattern at 8010 N. 
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Figure 19. Radiograph of two-ply panel at 9350 N ( 63% failure load). 



Figure 20. Brittle lacquer crack pattern at 9350 N showing 
transverse growth of severe damage region. 
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Figure 21. Radiograph of two-ply panel at 11,100 N (79% failure load). 







Figure 24. Brittle lacquer crack pattern at 12,500 N 
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Figure 25. Radiograph of two-ply panel at 12,800 N (93% failure load) showing 
additional fiber breaks at the right end of the notch. 



Figure 26. 


Brittle lacquer crack pattern at 12,800 N 
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Figure 27. 


Figure 28. Brittle lacquer crack pattern at 13,200 N 

showing large transverse progression of the 
severe damage region. 
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General Results and Discussion 


By examining the sequence of photographs and radio- 
graphs previously discussed, one can form a very clear pic- 
ture of the manner in which damage accumulates in a unidi- 
rectional laminate as the failure load is approached. 

Similar sequences could be presented for the entire range of 
specimens since the same testing procedure was followed for 
every specimen used in the study. The two-ply sequence 
presented, however, was quite typical of all of the tests 
and it was thus felt that the remainder of the results could 
be presented more effectively in the form of graphs which 
show results of all of the tests. Each of the questions 
raised in the introduction will be answered and discussed 
thoroughly in the remainder of this section. 

The question of thickness effect on the behavior of a 
unidirectional laminate can now be answered by examining the 
results of each set of tests. A full range of specimen 
thicknesses was tested for each panel width and initial 
notch width. It was observed that the severe matrix damage 
region as well as the overall appearance of the brittle 
coating at any given remote strain was very similar for the 
single-ply, the two-ply, the four-ply, and the eight-ply. 

The maximum length of this severe longitudinal matrix damage 
region (at failure) versus the number of plies is plotted 
for various numbers of crack lengths in Figure (29). Excel- 
lent agreement between the length of the observed damage 
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region and that predicted by the model exists as discussed 
in the previous section. This behavior is also discussed 
later in this section when a comparison is made between 
laminates having different matrix properties. 

It was also noted that the dimensions of the transverse 
damage region were very similar for a given initial notch 
width. The total number of broken fibers (per ply) ahead of 
the notch was approximately the same for the one-, two-, 
four-, and eight-ply laminates, as one can see in Figure 
(30). Thus, from these two graphs, is it clear that the 
longitudinal yielding and transverse damage regions depend 
very little on the number of plys in the unidirectional 
laminate for relatively thin laminates. 

Figure (31) shows the crack opening displacement (COD) 
for the set of specimens with an initial notch of nineteen 
broken fibers (or 3.17 mm). It is clear that the COD for 
each thickness is almost the same for a given percentage of 
the failure load and therefore the same remote strain. It 
can also be noted that the first fiber break occurs between 
50 and 55 percent of the ultimate failure load for each 
thickness tested. This phenomenon was observed in all of 
the tests regardless of notch width or specimen thickness. 
Figure (31) also indicates that the maximum COD at failure 
was almost identical for the four thicknesses tested. 

Figures (32) and (33) show the COD behavior for speci- 
mens with 53 and 89 initially broken fibers, respectively. 
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It is again clear that the specimen thickness has little 
effect on the COD for a given percentage of the failure 
load. The first fiber breaks occur at roughly the same 
remote strain level for a given notch width, as shown in the 
earlier graph. The appearance of the brittle lacquer coat- 
ing was observed to be the same for each thickness as in 
earlier tests. Also, the stress necessary to break the 
first fiber, the length of the severe damage region, and the 
notched strength of each laminate all are in close agreement 
with the quantities predicted by the mathematical model. 

This agreement will be discussed in detail later in this 
section. The similarity in behavior regardless of the spec- 
imen thickness (for up to eight plies) was observed for all 
initial notch widths, thus clearly indicating that a thick 
panel with unidirectional fibers in each ply can indeed be 
modeled as an equivalent monolayer. 

Another important objective of the present study was to 
determine the effect of different matrix properties on the 
fracture behavior of unidirectional laminates. Toward this 
end, certain sets of panels (as detailed in Table I) were 
heat-treated to bring their matrix properties close to the 
6061-T6 condition instead of the fully annealed 6061-0 con- 
dition as in the other panels. It was felt that a matrix 
with a different effective stiffness and yield strength 
would alter the fracture behavior since it is the matrix 
alone that transfers load to adjacent fibers in the lami- 
nate. 
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Figure (34) clearly shows that again the COD at a given 
percentage of the ultimate failure load is not strongly 
dependent on the number of plys in the laminate for the 
6061-T6 panels. As in the 6061-0 tests, the first fiber 
break occurs at 50 to 55 percent of the failure load. The 
COD at failure for the four thicknesses is almost the same, 
though it should be noted that this COD is much less than 
the COD for a similar panel with an annealed aluminum 
matrix. This change will be discussed fully, as it will 
help to determine the effect of matrix properties on the 
damage growth. 

Figure (35) is a remote stress-COD curve for various 
notch widths in a four-ply, 6061-0 matrix panel. As 
expected, specimens with short notches are stiffer than 
specimens with longer notches. Figure (36) shows the COD 
behavior for the same notch widths in a four-ply, 6061-T6 
matrix panel. It is clear that there is a drastic differ- 
ence in the notch opening under load when the properties of 
the matrix are changed. The stronger 6061-T6 matrix does 
not allow the fibers to displace nearly as much as the fully 
annealed 6061-0 matrix. 

It was suggested in [2] that the mathematical model 
employing the shear-lag assumption could predict accurately 
the load-COD curves for specimens of various notch widths if 
certain laminate material properties were known. Specifi- 
cally, these are the normalized matrix yield stress and the 
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effective shear stiffness of the fiber-matrix region which 
are calculated from an experimental load-COD curve. By 
using these properties, the load-COD behavior for any number 
of broken fibers could then be predicted accurately. Agree- 
ment was demonstrated in [2] for the narrow specimens and 
short notch lengths presented by Awerbuch and Hahn [6]. 

Using the much more extensive results of the present study 
it can be shown that, until large transverse damage is pres- 
ent, the load-COD curves for longer notch widths can also be 
matched. 

In Figure (37) a typical load-COD curve for a 4-ply 
unidirectional boron/aluminum composite panel is given for 
specimens having an overall width of 73.1 mm and 89 ini- 
tially broken fibers. Both a specimen having a 6061-0 ann- 
ealed matrix and a 6061-T6 matrix are shown, along with the 
load-COD curves predicted by the model. Using the values 
for the yield stress, t , and the effective shear stiffness, 
G^/h, given in [2], the load-COD curve for the 6061-0 speci- 
men can be closely matched at lower load levels even for 
longer notches than in [6]. The agreement is not as good 
for large notch extension as discussed later in this sec- 
tion. It must now be determined if appropriate changes in 
the two properties mentioned above will allow the model to 
predict accurately the load-COD curve for the 6061-T6 speci- 
men. 


76 



The load-COD curve for the 6061-T6 panel can be matched 

by increasing the value of the matrix yield stress 3.5 times 

6 2 

to 308 X 10 N/m and the value of the effective shear 

12 3 

stiffness 1.5 times to 172.5 X 10 N/m . These changes not 
only seem quite reasonable but the yield stress change 
agrees well with standard aluminum properties. The change 
in the value for Gj^/h is attributed more to a change in the 
shear transfer length, h, than to a change in G^ which 
should remain' relatively constant with any heat- treatment . 
The yield stress for 6061-T6 aluminum, on the other hand, is 
roughly 3-4 times the yield stress of 6061-0 aluminum. 

The results given in Figures (35), (36), and (37) for a 

given number of broken fibers are seen to be the same for 
lower loads, (i. e., below that required for matrix yielding 
at the crack tip) . In this low load range the experimental 
COD results are also predicted accurately by the orthotropic 
elasticity solution for a crack as given by Lekhnitskii 
[18 ] • For the elasticity solution the results are linear 
and the shear modulus, G^, was taken to be the same for both 
matrix materials. The 6061-0 matrix follows this linear 
representation for only the first few percent of the total 
load while the 6061-T6 matrix remains linear over about ten 
to fifteen percent of the total load. As the material model 
used in the analysis [2] for the matrix was taken to be . 
elastic-perf ectly plastic it is reasonable to expect that 
the model would fit the experimental results better for the 
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6061-T6 than for the 6061-0. This is indeed indicated by 
figure (37). The fact that the 6061-0 calculated results 
have a different initial slope than the experimental results 
also follows from this bi-linear fit to the' highly non-lin- 
ear behavior of the annealed matrix. 

The major deviation from the predicted curves occurs 
after several fibers have broken and the notch has begun to 
extend a significant amount. The extent of transverse notch 
extension for various notch widths and matrix tempers is 
shown in Table (IV) below. 


Table IV» Extent of stable notch extension at failure for 
selected initial notch widths in B/Al panels. 



2a=3.18 mm 

2a=9»52 mm 

2a=l 5*9 tnm 

Matrix 

W=25.4 mm 

W= 50.8 mm 

W=?3.1 mm 


Constrained Fibers 3 

4 

5 

6 061-0 

Total Number of o 

Broken Fibers J 

8 

14 


Constrained Fibers 2 

3 

3 

6061-T6 

Total Number of 2 

5 

7 


Broken Fibers 

( 


It was observed in all cases that the notches in the 
specimens having annealed matrices ran much further before 
fracture than the :6061-T6, Matrix yielding and broken fibers 
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accentuate or increase the notch opening. Thus, the COD for 
the annealed matrix specimens will be considerably greater 
than the COD for the T6 matrix at a load level corresponding 
to significant damage. It should be noted that the number 
of constrained fibers around which the aluminum matrix is 
still intact is almost constant for all initial notch widths 
for a given matrix temper. This number of constrained 
fibers also agrees with the number required by the mathe- 
matical model to fit the strength results. 

It was mentioned in the previous section that the 
appearance of the brittle coating indicated that there may 
be more than one yield region near the end of the notch. As 
each successive fiber break occurs, the newly broken fiber 
can move with respect to the unbroken adjacent fiber. Thus 
a new region of significant displacement and matrix yielding 
will form, allowing the COD to increase more than if the 
matrix remains elastic. At the present time the models of 
12 ] can incorporate only one yield region between the last 
broken fiber and the first unbroken fiber. Thus, if more 
than one yield region forms, the model cannot predict accu- 
rately the displacements of fibers after the notch extension 
has begun. 

Future plans include the modification of the mathe- 
matical model to include a new yield region as each successive 
fiber is broken during a loading cycle. However, substan- 
tial information can be obtained by summarizing the present 
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experimental data. Figure (38) shows the crack opening dis- 
placement at failure as a function of the number of ini- 
tially broken fibers. Clearly the crack opening displace- 
ment for the 6061-T6 panels is substantially less than that 
of the fully annealed 6061-0 panels. In fact, the differ- 
ence becomes larger as the number of initially broken fibers 
is increased. This phenomenon is caused by the relatively 
large notch extension of the fully annealed specimens. The 
stronger T6 matrix cannot absorb as much energy through 
yielding as the annealed matrix, thus the load is transfer- 
red to the next unbroken fiber sooner and the notch exten- 
sion becomes unstable, leading to earlier failure. 

The longitudinal damage length as described in an ear- 
lier section is plotted as a function of the number of bro- 
ken fibers in Figure (39). The predicted value for this 
damage length is obtained by calculating the yield length 
necessary to reach the breaking stress in the next unbroken 
fiber through the use of the mathematical model. This value 
should correspond to the severe damage region on the experi- 
mental specimens since there should be no further yielding 
in the region between fibers once the next fiber has broken 
and is free to displace. It is clear that the predicted 
trend of damage growth from the model is a very good indica- 
tor of the actual damage present in a unidirectional lami- 
nate. The stronger T6 matrix substantially reduces the 
yielding in the matrix and thus reduces the amount of damage 
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tolerance possessed by a boron/aluminum laminate. The 
result of increasing the yield strength of the matrix mate- 
rial used in a unidirectional laminate is thus to make the 
laminate more notch sensitive and actually weaker than a 
laminate having a more ductile, energy-absorbing matrix. 

The notched strengths of the test panels and the reduc- 
tion in strength for the 6061-T6 matrix can be studied in 
Table (V) and in Figure (40). It should again be noted that 
the effect of the change in matrix properties becomes much 
more evident as the notch becomes larger; that is, the num- 
ber of initially broken fibers is increased. As in [2], the 
strength curve for the fully- annealed 6061-0 matrix speci- 
mens can be predicted by the mathematical model using 6 con- 
strained fibers ahead of the notch. It should be remembered 
that in the damage photograph sequence shown earlier there 
were approximately 6 constrained fibers ahead of the notch 
tip even for large notch extension. The strength curve for 
the 6061-T6 panels can be approximated by using 3 con- 
strained fibers ahead of the notch. This change has physi- 
cal significance; the stable notch extension is larger for 
the fully annealed specimens than for the heat-treated spec- 
imens. Thus as discussed earlier, the heat-treating of the 
matrix actually reduces the damage tolerance capabilities of 
the laminate and makes the unidirectional laminate weaker 
than one with a fully annealed matrix. It is noted that the 
notched strengths found in this study for the 6061-T6 matrix 
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agree well with notched strengths found by Awerbuch and Hahn 
in [ 19 ]. 

One should be careful in extending this reasoning to 
even stronger matrix materials in which the transition is 
made from longitudinal yielding to longitudinal splitting. 
For example, an epoxy matrix is one in which little yielding 
is exhibited but longitudinal splits occur at the notch tips 
under low load. This phenomenon is discussed in [ 2 ] and has 
been observed by many investigators. It is important to 
realize that once this small split forms, it will quickly 
grow to the machine grips in an experimental test (or to 
infinity as predicted by the model). The stress concentra- 
tions are then removed and the notched laminate will fail at 
the net section failure load. Thus the conclusion that a 
stiffer or stronger matrix will weaken the notched laminate 
is valid only for matrices that are not brittle enough to 
form longitudinal splits between fibers at the end of the 
notch. Some small splitting was observed in the present 
study in 6061-T6 laminates with very long notches, possibly 
indicating that the T6 matrix is close to the transition 
from matrix yielding to splitting. 
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SEVERE DAMAGE LENGTH (cm) 



Figure 29. Variation of the severe damage length with the 
number of plies for selected numbers of broken 
fibers . 



NUMBER OF BROKEN FIBERS AHEAD OF NOTCH 



Figure 30. Variation of the transverse notch extension 

with the number of plies for selected numbers 
of broken fibers. 



PERCENT FAILURE LOAD 



CRACK OPENING DISPLACEMENT (yjn) 

Figure 31. Percent failure load versus crack opening displacement for 19 broken fibers 
per ply; comparison of various numbers of plies for a 6061-0 A1 matrix. 
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displacement for 35 broken fibers 
■ plies for a 6061-0 A1 matrix. 




Figure 33. Percent failure load versus crack opening displacement for 89 broken fibers 
per ply; comparison of various numbers of plies for a 6061-0 A1 matrix. 



PERCENT FAILURE 



Figure 34. 


Percent failure load v 
per ply; comparison of 


NOTCH WIDTH (2a) = 15.9 mm 
89 BROKEN FIBERS PER PLY 


FIRST FIBER 
BREAK 
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A.CK OPENING DISPLACEMENT (pi) 

•sus crack opening displacement for 89 broken fibers 
r arious numbers of plies for a 6061-T6 A1 matrix. 



REMOTE STRESS (MPa) 




REMOTE STRESS (MPa) 



CRACK OPENING DISPLACEMENT (urn) 


Figure 36. Remote stress versus crack opening displacement 
for selected numbers of broken fibers in panels 
having a 6061 -T6 aluminum matrix. 


LOAD (kN) 



CRACK OPENING DISPLACEMENT (vim) 


Figure 37. Load versus crack opening displacement for 89 
broken fibers in a 4-ply panel; comparison of 
the analytical and experimental curves . 
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COD AT FAILURE (p) 



Figure 38. Crack opening displacement at failure as a function of the number of broken 
fibers; comparison of the different matrix tempers. 


DAMAGE LENGTH (cm) 



Figure 39. Severe damage length at first fiber break as a function of the number of 
broken fibers; comparison of the different matrix tempers. 
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0 20 40 60 80 100 120 

NUMBER OF BROKEN FIBERS 

Figure 40. Strength curve: Applied stress as a function of the number of broken 

fibers; comparison of the different matrix tempers. 


CHAPTER V 


CONCLUSIONS 

Results of this study indicate that the mathematical 
models presented in [2] based on the classical shear-lag 
assumption predict the displacements and stresses in a uni- 
directional fibrous composite laminate with very acceptable 
accuracy for a wide range of notch lengths, laminate widths, 
and laminate thicknesses. It is also shown that the crack 
opening displacement versus load curve predicted by the 
model is very close to the actual behavior, particularly in 
the early part of the load cycle before a large number of 
fibers have broken. The extent of the stable notch exten- 
sion, that is, the number of broken but constrained fibers, 
predicted by the model is found to accurately represent that 
observed in all laminates tested. Also, the dimensions of 
the longitudinal damage region observed in the laboratory 
are found to agree well with those calculated by the shear- 
lag model, as were the notched strengths of the test speci- 
mens . 

It is demonstrated that the thickness (number of plies) 
has very little effect on the damage growth and fracture 
behavior of a unidirectional composite laminate. Longitudi- 
nal damage regions and crack opening displacements measured 
for single-ply, two-ply, four-ply, and eight-ply specimens 
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are found to be almost identical for a given value of remote 
strain. The stable notch extension in the form of broken 
fibers and matrix yielding ahead of the notch is found to be 
present in the same quantities in specimens having various 
numbers of plies. 

Finally, it is shown that heat-treating the matrix to 
increase the yield strength actually reduces the notched 
strength of the laminate. The stronger matrix transfers 
load to the unbroken fibers at a lower load than the fully 
annealed matrix, with a substantial reduction in the amount 
of matrix yielding. The results of these tests indicate 
that to increase the capability of a laminate to tolerate a 
defect during a load cycle, the matrix must be able to yield 
and thereby absorb some energy that would otherwise serve to 
propagate the notch and cause failure of the laminate. 
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